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THE BLACK HOLE MASSES AND HOST GALAXIES OF BL LAC OBJECTS 
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ABSTRACT 

We have measured the central stellar velocity dispersion in the host galaxies of 1 1 BL Lac objects with redshifts 
z < 0.125. The range of velocity dispersions, ~ 170-370 km s" 1 , is similar to that of nearby radio galaxies. 
Using the correlation between stellar velocity dispersion and black hole mass defined for nearby galaxies, we 
derive estimates of the black hole masses in the range 10 7,9 - 10 92 M Q . We do not find any significant difference 
between the black hole masses in high-frequency-peaked and low-frequency -peaked BL Lac objects. Combining 
the velocity dispersions with previously measured host galaxy structural parameters, we find that the host galaxies 
lie on the fundamental plane of elliptical galaxies. This supports the conclusions of imaging studies that the 
majority of BL Lac hosts are normal giant ellipticals. 

Subject headings: BL Lacertae objects: general — galaxies: active — galaxies: elliptical and lenticular — 
galaxies: kinematics and dynamics — galaxies: nuclei 



1. INTRODUCTION 

One of the goals of AGN research is to develop a unified 
framework in which the diversity of the AGN family might be 
understood in terms of variations in a few fundamental parame- 
ters such as the black hole mass, the ratio of the accretion rate to 
the Eddington accretion rate, and the orientation relative to our 
line of sight (e.g., Laor 2000; Boroson 2002). Determining the 
masses of the central black holes in different classes of AGNs is 
an important step toward this goal. An key advance in this area 
was the recent discovery of the Mbh - cr relation, a tight corre- 
lation between black hole mass and stellar velocity dispersion 
in the host galaxy bulge (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000a). The Mbh-ct relation has provided an important 
consistency check for black hole masses determined from re- 
verberation mapping of Seyfert nuclei (Nelson 2000; Gebhardt 
et al. 2000b; Ferrarese et al. 2001). The tightness of the correla- 
tion over a wide range of host galaxy types, both elliptical and 
spiral, makes it tremendously useful as a means to determine 
the black hole masses in AGNs. While stellar-dynamical or 
gas-dynamical black hole mass measurements are particularly 
difficult for AGNs, stellar velocity dispersions can be measured 
relatively easily in some classes of active galaxies. Further- 
more, black hole mass estimates based on velocity dispersions 
are expected to be much more accurate than estimates derived 
from the loose correlation between Mbh and bulge luminosity 
(Magorrian et al. 1998; Kormendy & Gebhart 2001). 

We recently began a program to measure the stellar velocity 
dispersions in the host galaxies of low-redshift BL Lac objects 
in order to apply the Mbh _ c relation and determine their black 
hole masses. The first measurement of a velocity dispersion in 
a BL Lac object, for the TeV 7-ray source Markarian 501, in- 
dicated a likely black hole mass of (0.9-3.4) x 10 9 M Q (Barth, 
Ho, & Sargent 2002a, Paper I). In this paper, we report on mea- 
surements for a sample of 1 1 objects and discuss the results 
in the context of unification models for radio-loud AGNs. We 
also compare our results with measurements of stellar veloc- 
ity dispersions for BL Lac objects reported recently by Falomo, 
Kotilainen, & Treves (2002, hereinafter FKT). 



2. OBSERVATIONS AND REDUCTIONS 

The spectra were obtained at the Hale 5m telescope at Palo- 
mar Observatory, and at the 2.5m du Pont and 6.5m Baade 
(Magellan 1) telescopes at Las Campanas Observatory. We 
used the Double Spectrograph at Palomar (Oke & Gunn 1982), 
the ModSpec Spectrograph at the du Pont telescope, and the 
Boiler & Chivens Spectrograph at Magellan. Details of the in- 
strumental setups are given in Table 1; Table 2 describes the 
exposures for each object. Conditions were mediocre during 
the Palomar nights, with typical seeing of 1 ."5 to 2."0 and thick 
clouds at times, while the Las Campanas nights were clear. 
The gratings and grating angles were chosen to cover the most 
useful spectral regions for measuring velocity dispersions: the 
Mg lb region in the blue, and the Ca II near-infrared triplet in the 
red. The values for the seeing listed in Table 2 were measured 
from observations of template and standard stars, but these were 
mostly observed during twilight and the seeing during the BL 
Lac observations may have been somewhat different. For each 
galaxy, the slit was oriented close to the parallactic angle at the 
midpoint of the observation, and the total exposure time was 
divided into individual exposures of typically 1800 s. 

Bias subtraction, flat-fielding, spectral extractions, and wave- 
length calibration were done using standard routines in IRAF, 
and the flux calibration and subsequent measurements were per- 
formed in IDL. Spectra were extracted using the optimal algo- 
rithm of Home (1986) as implemented in IRAF; the extraction 
widths are listed in Table 2. The strongest telluric absorption 
bands were removed by division by the normalized spectrum 
of a nearly featureless star following methods similar to those 
of Wade & Home (1988). After wavelength calibration, a fi- 
nal shift to the wavelength scale was applied to each exposure 
based on the wavelengths of night sky emission lines. This 
step ensured that the stellar absorption features were not arti- 
ficially broadened by summing exposures with slightly offset 
wavelength scales. 

During each observing run, several stars were observed for 
use as velocity templates. The template stars were chosen to be 
primarily giants of spectral type between G8 and K5, but we 
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also observed some subgiants and dwarfs and also a few F and 
early-G stars. In the following analysis, we consider only the 
G and K giant templates, since they give the best match to the 
dominant stellar populations in the blazar host galaxies. 

In addition to the sample listed in Table 1, we also observed 
the object BL Lac at Palomar, but were unable to measure its ve- 
locity dispersion. The stellar absorption lines were only barely 
detected in a total exposure of 7200 s due to overwhelming 
AGN contamination. Also, BL Lac has a Galactic extinction 
of Ag = 1.42 mag (Schlegel, Finkbeiner, & Davis 1998), and 
the dominant features over much of its optical spectrum are the 
Galactic diffuse interstellar bands (Vermeulen et al. 1995). 

3. MEASUREMENT OF VELOCITY DISPERSION 

Measurement of the velocity dispersion of a BL Lac host 
galaxy is straightforward in principle, but there are some com- 
plications arising from the dilution of the host galaxy starlight 
by nonthermal emission. For the observations presented here, 
the detected flux within the spectroscopic aperture is ~ 50 - 
80% nonthermal continuum emission. As a result, a BL Lac ob- 
ject spectrum yields a more uncertain velocity dispersion than a 
normal galaxy spectrum of the same signal-to-noise (S/N) ratio. 
Furthermore, Galactic or telluric absorption features, if not cor- 
rected for, will have a relatively larger impact on the velocity 
dispersion measurements for the BL Lac objects. 

We measured the velocity dispersions of the BL Lac hosts 
by fitting broadened stellar spectra to the object spectra, fol- 
lowing a method similar to those described by van der Marel 
(1994), Rix et al. (1995), and Kelson et al. (2000). The object 
and template spectra are first rebinned to a wavelength scale 
that is linear in log A, and the object spectrum is transformed to 
zero redshift. The template spectrum is broadened by convo- 
lution with a Gaussian, and adjusted to match the galaxy spec- 
trum by adding a linear featureless continuum and multiplying 
by a low-order polynomial. The fit includes six free parameters: 
the systemic velocity of the galaxy, the width of the Gaussian 
broadening function, the normalization and slope of the feature- 
less continuum, and the three coefficients of a quadratic polyno- 
mial. Multiplication by the polynomial allows the fit to account 
for any reddening or other differences in overall spectral shape 
between the galaxy and the template, without introducing any 
high-frequency structure that might affect the dispersion mea- 
surement. The fit is optimized using a downhill simplex algo- 
rithm (Press et al. 1988). For each template, the fit is repeated 
several times with a range of initial parameter values to ensure 
that the global minimum of \ 2 is found. The la uncertainty 
for each template is determined by finding how far the velocity 
dispersion must be varied from its best-fit value (with all other 
parameters allowed to float freely) in order to increase x 2 by 
unity relative to its minimum value. 

The wavelength range for the fits was selected to include the 
regions most sensitive to the stellar velocity dispersion. For 
each object, we attempted to use the largest wavelength range 
over which a good fit could be achieved, excluding emission 
lines or regions that consistently gave poor fits with a range of 
template stars. In general, velocity dispersions measured from 
the Ca triplet lines are preferred, if they can be observed, be- 
cause the Ca triplet is less sensitive to template mismatch than 
the blue spectral region (Dressier 1984). Also, the Ca triplet 
lines are intrinsically strong and well isolated from neighboring 
lines, and can be detected even in highly diluted spectra. For the 
lowest-redshift objects (Mrk 421, Mrk 501, and 2201+044), we 



used the Ca triplet lines to measure a. 

At recession velocities beyond ~ 9000 km s" 1 , the Ca II 
A8662 line is redshifted into deep telluric water vapor absorp- 
tion bands that begin at about 8925 A. Reliable velocity dis- 
persions can still be measured using the Ca lines at 8498 and 
8542 A out to cz ~ 12000 km s" 1 . For the objects having larger 
redshifts, we used the region around Mg lb to measure the ve- 
locity dispersions. Best results were obtained for the region 
5205 - 5450 A, which contains reasonably strong features due 
to Ca and Fe. We also included a wavelength range just blue- 
ward of Mg lb in the fit, typically 5030-5150 A, for four ob- 
jects having sufficient S/N in this region. The wavelength re- 
gion was adjusted for some objects to exclude Galactic absorp- 
tion features or bad pixels, or to include a wider wavelength 
range when a good fit could be obtained. 

It proved impossible to fit both the 5205 - 5450 A region 
(which contains the Fe 5270 and Fe 5335 blends) and the Mg lb 
line simultaneously. This is a consequence of the correlation 
between the [Mg/Fe] ratio and velocity dispersion in elliptical 
galaxies (Worthey, Faber, & Gonzalez 1992; Trageret al. 1998; 
Kuntschner et al. 2001). Since the Mg lines increase in strength 
more steeply than do the Fe lines as a function of er, it is difficult 
to find local template stars that will provide an exact match to 
the [Mg/Fe] ratio of any individual giant elliptical galaxy. This 
source of template mismatch is a particular problem for direct 
fitting measurements, since the calculation of \ 2 is sensitive to 
both the width and the depth of absorption lines. Thus, we ex- 
cluded the Mg lb line from the fits, preferring to fit the larger 
surrounding region. Barth et al. (2002b) discuss the [Mg/Fe] 
template matching problem in further detail, using spectra of 
nearby galaxies obtained during the same observing runs. The 
best fits for each object are shown in Figure 1. 

The dominant contribution to the measurement uncertainty in 
most cases is template mismatch, since the scatter in a measure- 
ments for different templates is generally larger than the fitting 
uncertainty for a single template. The final velocity dispersions 
listed in Table 2 are the weighted mean of values measured 
with several templates of type G8-K5. To account for tem- 
plate matching errors, the measurement error for each galaxy 
is the sum in quadrature of the fitting error for the best-fitting 
template and the standard deviation of the measurements for all 
templates. 

There is one additional source of systematic uncertainty that 
contributes significantly to the error budget. As discussed by 
Barth et al. (2002b), tests on high S/N spectra of nearby galax- 
ies demonstrate that our fitting results for the Mg lb region are 
somewhat sensitive to the order of the polynomial that is used 
to match the continuum shape of the template to the galaxy. 
For polynomial orders of 2, 3, or 4, we find variations in the 
derived value of a at the 5% level. To account for this system- 
atic uncertainty in the model fits, we add 0.05cr in quadrature 
to the uncertainties on each measurement of a derived from the 
blue spectral region. The Ca triplet measurements are much less 
sensitive to the polynomial order, with variations at the level of 
^ 1%. 

Another potential source of systematic error is the variation 
of the instrumental linewidth across the detector, since a given 
spectral feature will fall on a different portion of the CCD 
for different redshifts. This would have the largest effect on 
the Palomar blue side measurements with the 600 line grating, 
which was the setting with the lowest spectral resolution, and 
also with the largest focus variations across the detector. From 
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examination of comparison lamp spectra, we estimate that the 
maximum variation in instrumental dispersion over the region 
used for the measurements is <~ 15 km s" 1 in this setting. Since 
the instrumental dispersion contributes in quadrature to the total 
observed linewidth, these variations will have a very small im- 
pact on the final measurements. For the higher-dispersion set- 
tings the impact of focus variations on the measurements will 
be negligible. 

3.1. Tests of the Measurement Routine 

At Palomar, we also observed six velocity dispersion stan- 
dard galaxies from the catalog of McElroy (1995), with a rang- 
ing over ~ 100-300 km s" 1 , as a check on the fitting routine. 
The results are listed in Table 3. The blue and red side measure- 
ments for each galaxy agree within their la uncertainty ranges. 
We find good agreement between our measurements and the av- 
erage velocity dispersion data compiled by McElroy for four of 
the six galaxies, while our dispersions for NGC 4579 and NGC 
4736 are significantly lower than the McElroy average values. 
However, in both cases our measurements for these objects do 
fall in the range of previously measured values as listed in the 
Hypercat database, 4 and the close agreement between the blue 
and red measurements suggests that our values of a do give 
an accurate measure of the velocity dispersion within the ob- 
served aperture. We believe that the differences between our 
results and the mean catalog results for these galaxies are due 
to individual, discrepant measurements that are included in the 
catalog average velocity dispersions. A full description of the 
nearby galaxy measurements is given by Barth et al. (2002b). 

Our observations have a wide range in S/N, and one possi- 
ble concern is that at low S/N the fitting routine could be bi- 
ased in a systematic way. To test the accuracy of the results, 
we performed Monte Carlo simulations to determine how well 
we could recover the correct velocity dispersion as a function 
of S/N. Starting with the Palomar blue side spectrum of NGC 
2841, which was observed with the 1200 line grating, we added 
a featureless f\ oc A" 1 continuum normalized to contribute 2/3 
of the total flux at 5500 A. Then, Poisson noise was added to 
the spectrum to achieve a desired S/N ratio. For each value of 
the S/N ratio (ranging from 500 to 50 per pixel), we created 100 
realizations of the diluted, noise-added spectrum and measured 
their velocity dispersions using the same measurement routine 
and fitting over the wavelength ranges 5050-5150 and 5205- 
5450 A. Ten template stars were used for each measurement. 

Measurement of the original spectrum gives a = 229 ± 5 km 
s" 1 , in good agreement with previous results from the literature 
(McElroy 1995). The results of the Monte Carlo calculations 
are displayed in Figure 2. At each S/N ratio, we computed the 
weighted average and the standard deviation of the velocity dis- 
persion measurements for the 100 realizations of the spectrum. 
Over the full range of S/N used in the calculations, the mean 
velocity dispersion is within ±3 km s" 1 of the expected result 
(229 km s" 1 ), demonstrating that there is no systematic offset 
as a function of S/N. As expected, the scatter in the measure- 
ments rises at low S/N, but for S/N > 100 per pixel the standard 
deviation is less than 15 km s" 1 . All but one of our observations 
have S/N > 100 per pixel, so the expected measurement error 
is not large. The actual expected error for any measurement is 
a function of the S/N, dilution, spectral resolution, and velocity 
dispersion, and we do not attempt to calculate it in more detail, 
but the simulations demonstrate that the template fitting routine 

4 The Hypercat database can be found at http://www-obs.univ-lyonl.fr/hypercat/ 



should work well for our data. 

3.2. Notes on Individual Objects 

Mrk 421: The Ca triplet lines are easily detected and we 
found a good fit to the A8542 A profile for a = 216 km s _1 . The 
A8662 line was not used in the fit as the telluric H2O absorp- 
tion beyond 8900 A could not be adequately removed from the 
spectrum by standard star division. The blue side spectrum was 
more strongly diluted by nonthermal light than the other galax- 
ies in the sample, and it proved difficult to obtain an adequate 
fit. For five template stars, the blue side data yields a = 2 1 3 ± 47 
km s" 1 . Thus, we use the red side measurement as our final re- 
sult. 

Mrk 501: The Ca II triplet region gives a = 372 ± 18 km 
s" 1 , where the uncertainty range includes the scatter among 
24 template stars ranging between F3 and K5 and including 
giants, subgiants, and dwarfs (Paper I). In the blue spectral 
range, the best-fitting template (HD 125560, a K3III star) gives 
a = 386 ± 9 km s" 1 over 5200-5600 A, and the standard devia- 
tion of results with all templates is 8 1 km s" 1 . However, if we 
restrict the fitting to only G and K giants, the scatter is reduced 
to ±11 kms" 1 . 

During the night that Mrk 501 was observed, the blue side 
flatfield structure shifted on the detector along the dispersion 
direction, presumably due to mechanical flexure in the spec- 
trograph. We observed dome and internal flats in the afternoon. 
Objects observed during the first half of the night flattened well, 
but in objects observed during the second half of the night (in- 
cluding Mrk 501) there were some residual bumps in the flat- 
tened images at the ~ 10% level at A < 5000 A. Although this 
did not visibly affect the region used to measure er, and the blue 
and red side measurements are in agreement, we consider the 
red side measurement to be more reliable and we use it as our 
best estimate of a. 

A Hubble Space Telescope (HST) image of Mrk 501 was 
taken as part of the snapshot survey of Scarpa et al. (2000) 
and Urry et al. (2000), but it was not included in the analysis 
of their HST sample. To measure the structural parameters of 
the host galaxy, we downloaded the images from the HST data 
archive. There are 5 separate exposures with exposure times of 
2 s, 2 x 30 s, and 2 x 120 s, all with the F702W filter (similar 
to ground-based R). The nucleus is saturated in the 30 s and 
120 s exposures. We combined the 120 s exposures, eliminat- 
ing cosmic -ray hits, and used the GALFIT program (Peng et al. 
2002) to perform a two-dimensional fit to the host galaxy pro- 
file. The model used in the fit was a de Vaucouleurs (1948) 
profile combined with an unresolved nuclear point source, con- 
volved with a synthetic point-spread function (PSF) generated 
by the Tiny Tim program (Krist & Hook 2001). Saturated pix- 
els at the nucleus were masked out in the fits to the longer 
exposure. For consistency with the Urry et al. (2000) sam- 
ple, we converted the F702W filter magnitude to the Cousins 
R band with the SYNPHOT package in IRAF, using a redshifted 
elliptical galaxy spectrum to calculate the correction. The total 
host galaxy magnitude is m« = 13.07 mag, in good agreement 
with the ground-based measurement of Wurtz, Stocke, & Yee 
(1996). The fit gives an effective radius of r e = 4"4, substan- 
tially smaller than previous ground-based measurements (Nils- 
son et al. 1999); the comparison between ground-based and 
HST measurements of r e is discussed further in §5. 
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0521-365 and 0548-322: During the Magellan run, only one 
K-giant stellar template was observed (HR 1629, K4III). To al- 
low for template mismatch error, we add 20 km s" 1 in quadra- 
ture to the uncertainty in the velocity dispersion measurement 
from this template. 

0706+591: At z = 0. 125, this is the most distant object in the 
sample. Galactic Na D absorption falls in the region used to 
measure a and is excluded from the fit. The broad H2O telluric 
absorption band between 5870 and 6000 A also falls in this re- 
gion. We could not entirely exclude this band from the fit as it 
covers the most useful portion of the spectrum; inclusion of part 
of the telluric band in the fitting region is probably responsible 
for the relatively poor fit in the neighborhood of the Mg lb line. 
The measurement yields a = 216 ± 23 km s" 1 . As a check on 
this result, we also performed a fit over the region 4920-5150 
A (excluding a 20 A region around [O III] A5007 although this 
line was not clearly detected), and found a = 204 ± 27 km s" 1 . 

2201+044: The core of the Ca II A8542 absorption line falls 
on a night sky emission line and the extracted spectrum has a 
sharp negative "spike" at this wavelength. We masked out a 
5 A-wide region surrounding this spike in the final fit, but the 
measured value of a changes by less than 10% if the entire line 
is included in the fit. Over the wavelength range used for the 
fit there are several extremely weak telluric absorption features 
with expected equivalent widths less than 0.02 A. They were 
not masked out in the fit, since their amplitudes are smaller than 
the noise level in the data. 

3.3. Comparison with the results ofFalomo et al. 

FKT reported velocity dispersion measurements for seven 
BL Lac objects, six of which are in our sample. Comparison 
of the results shows close agreement for some galaxies such 
as Mrk 421, but also some surprisingly large differences. The 
worst disagreement is found for Mrk 501 and 1 Zw 187. For 
these two galaxies, FKT find a c = 291 and 253 km s" 1 , respec- 
tively, while we obtain a = 372 and 171 km s _1 . FKT suggest 
that the different results for Mrk 501 could be the result of tel- 
luric absorption features in our red spectrum. However, this is 
not the case because we excluded these telluric lines from the 
fitting region, as described in Paper I. Different aperture sizes 
cannot be the cause of such large discrepancies either. FKT 
used a 1" slit and extraction widths of 3" -5", similar to the 
aperture sizes we used, and the resulting variations in a should 
be below the 10% level in most cases (Gebhardt et al. 2000a; 
Merritt & Ferrarese 2001). 

FKT used the Fourier quotient routine in IRAF which, in 
principle, should yield the same results on the same data as a di- 
rect fitting routine. However, their measurements were carried 
out over a much wider wavelength range containing different 
spectral features. They used two grating settings for their mea- 
surements: 4800-5800 A and 5700-8000 A (setups "A" and 
"B", respectively). These ranges include Ha and H/3, which 
are among the strongest absorption features in the spectrum but 
which also appear in emission in some objects. They also con- 
tain other possible emission features including [O III] A5007. 
(In our sample of 1 1 objects, [O III] emission is clearly vis- 
ible in five, and H/3 in three.) FKT removed emission lines 
from their spectral by connecting a linear continuum across the 
wavelength range spanned by the emission lines. This is a very 
poor approximation to the true starlight spectrum of the galaxy, 
particularly for the Balmer lines, and may adversely effect the 

5 See Fig. 16 of Matheson et al. (2000) for a plot of the telluric absorption spectrum 



results. It would have been preferable to replace the emission 
lines with the spectra of appropriately broadened template stars 
(van der Marel & Franx 1993) or with spectra of galaxies hav- 
ing similar velocity dispersion and line strength (Dalle Ore et al. 
1991), or instead to perform the fit over a wavelength range that 
is devoid of emission lines. Even when Ha and H/3 emission 
are not clearly visible in the spectra, they could still partially 
fill in the Balmer absorption lines and affect the measurements. 

The strongest absorption feature in the wavelength range 
used by FKT is the Na D doublet at 5890, 5896 A, but this 
could have a contribution from interstellar absorption in the 
host galaxies. There is Galactic interstellar Na D absorption 
in some objects as well; our spectra reveal strong Galactic Na 
D absorption in Mrk 421, and weak absorption in 0706+591. 
With the Fourier quotient routine there is no entirely satisfac- 
tory way to exclude such a feature from the measurement. Since 
there are few other strong stellar features within a few hundred 
A of Na D, this spectral region is considered to be one of the 
least useful portions of the optical spectrum for measuring ve- 
locity dispersions (Dressier 1984). The weak telluric O2 and 
H2O absorption bands between 5800 and 6600 A might also 
have a small effect on the measurements, especially in objects 
with a very high degree of dilution by nonthermal light. 5 

Whenever possible, our measurements were performed on 
regions selected to have good sensitivity to velocity dispersion 
without having significant contamination from nebular emis- 
sion, interstellar absorption, or telluric absorption. A major ad- 
vantage of the direct fitting method over the Fourier quotient 
method for these objects is that it provides direct visual feed- 
back on the quality of the template match. This makes it possi- 
ble to mask out regions which are not fit well by any template 
star, as well as telluric absorption bands and emission lines. We 
note, also, that our spectra have considerably higher S/N than 
those obtained by FKT; this enables us to perform our velocity 
dispersion measurements on a smaller spectral region. In sum- 
mary, we believe that our velocity dispersion determinations are 
more accurate than those of FKT. 



4. BLACK HOLE MASSES 

To determine M BH from the velocity dispersion, we use the 
recent fit to the M BH - relation from Tremaine et al. (2002), 
which is given by 

log(M BH /M Q ) = (8.13 ±0.06)+(4.02±0.32) log ( ^g— ^ . 

This relation was derived using a carefully culled sample of the 
best-quality BH mass measurements, and using a fitting tech- 
nique that properly accounts for uncertainties both in M BH and 
a. 

The measurements used to calibrate the above relation were 
luminosity-weighted velocity dispersions a e measured in a slit 
aperture of length 2r e . However, our values of o do not ex- 
actly correspond to this aperture size, and our spectra generally 
do not have sufficient sensitivity to measure a outside the inner 
few arcseconds. Velocity dispersions can be corrected to a stan- 
dard aperture size using the relations given by J0rgensen, Franx, 
& Kjaergaard (1995a), but these relations are not calibrated for 
long-slit apertures such as those used by Tremaine et al. (2002). 
Another difficulty is that the corrections depend on the effective 
radii, which are not known accurately for these objects. Rather 

over this wavelength range. 
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than apply an uncertain correction to the measured velocity dis- 
persions, we choose to use our measured dispersions directly to 
estimate M BH . As shown by Gebhardt et al. (2000a), the value 
of a measured from slit apertures of different length will differ 
from a e by less than 10% for typical, nearby elliptical galaxies. 

The results are listed in Table 2, with uncertainties in M SH 
determined by propagating the measurement uncertainty in a 
as well as the uncertainties in the coefficients of the M BH - o 
relation from Tremaine et al. (2002). We do not include the ad- 
ditional possible error due to intrinsic scatter in the Mbh - o re- 
lation, because the amount of scatter is not yet well determined. 
Tremaine et al. (2002) have determined that this intrinsic scat- 
ter is < 0.3 dex; it is hoped that future BH mass measurements 
with HST will provide a more definitive result. The masses 
range from 7 x 10 7 to 1.6 x 10 9 M Q , with a median value of 
1.8 x 10 8 M Q . 

We emphasize that these black hole masses are indirect esti- 
mates, not dynamical measurements of the mass. Nevertheless, 
there is no strong reason to believe that the Mbh _ f correla- 
tion would not apply to BL Lac objects. Imaging of BL Lacs 
with HST has demonstrated that the host galaxies are all appar- 
ently normal giant elliptical galaxies (Scarpa et al. 2000; Urry 
et al. 2000), in contrast to the disturbed morphologies seen in 
some QSO host galaxies (e.g., Bahcall et al. 1997). Thus, the 
stellar velocity dispersions in BL Lac objects should not be af- 
fected by recent mergers or other peculiarities. No nearby el- 
liptical galaxies have yet been found to deviate strongly from 
the Mbh _ cr relation, so it would be very surprising if the BL 
Lac hosts did not follow the correlation. Furthermore, the up- 
per end of the M B h - & relation is largely calibrated with FR I 
radio galaxies, including M87 and NGC 4261, which presum- 
ably would appear as BL Lac objects if viewed along the jet, 
and there is no evidence for any offset in the M B h - cr corre- 
lation of radio galaxies in comparison with inactive ellipticals. 
Still, it would be a useful consistency check to perform stellar- 
dynamical measurements of M B h in a few of the lowest-redshift 
BL Lac objects to confirm that these assumptions are indeed 
valid. 

5. DISCUSSION 

According to unified models for radio-loud AGNs, BL Lac 
objects are FR I radio galaxies viewed along the jet axis (e.g., 
Urry & Padovani 1995; Urry et al. 2000). If radio galaxies and 
BL Lac objects were derived from identical parent populations, 
then the distribution of black hole masses in BL Lac objects 
should be the same as that in nearby FR I radio galaxies. How- 
ever, different selection effects for BL Lacs and radio galaxies 
might be expected to bias the BL Lac objects toward less mas- 
sive and less luminous host galaxies. For example, Urry et al. 
(2000) point out that radio-loud AGNs with jets beamed in our 
direction might not be identified as BL Lac objects if they occur 
in very bright host galaxies so that their stellar continua are not 
diluted enough by nonthermal emission to match the standard 
defining criteria for classification as BL Lac objects. Thus, the 
naive expectation that BL Lac objects and radio galaxies should 
have identical distributions of black hole masses or other host 
galaxy properties is likely to be an oversimplification. Such dif- 
ferences do not violate the basic premise of unified models, as 
long as objects classified as BL Lacs would still appear to be 
normal radio galaxies when viewed from an oblique angle. 

As a comparison sample of nearby radio galaxies, we use the 
compilation by Bettoni et al. (2001, 2002). While this is not a 



complete sample in any sense, it is still the largest sample of 
nearby radio galaxies for which velocity dispersions, effective 
radii, and surface brightnesses are available. Figure 3 shows the 
distribution of velocity dispersions for the BL Lac objects and 
the radio galaxy sample of Bettoni et al. (2001). As demon- 
strated previously by FKT, the radio galaxies and BL Lac ob- 
jects occupy a similar range of values of cr, and therefore M BH . 
The two samples have nearly the same mean velocity disper- 
sion, 247 km s" 1 for the BL Lac objects and 256 km s" 1 for the 
radio galaxies. A Kolmogorov-Smirnov test (Press et al. 1988) 
on the two distributions shows that the hypothesis of identical 
parent populations for the two distributions can only be rejected 
at the 9% confidence level. Thus, the two distributions appear 
to be indistinguishable within the limited statistics, but since 
sample incompleteness and selection effects are not accounted 
for in any way, we caution that this comparison does not lead 
to any firm conclusions regarding the nature of the two popula- 
tions. 

Our sample consists of three low-frequency-peaked (LBL) 
and eight high-frequency-peaked (HBL) objects. Figure 3 illus- 
trates that, in agreement with FKT, there is no large systematic 
offset between the black hole masses in LBL and HBL objects. 
The LBL and HBL objects in our sample have <a> = 238 and 
235 km s" 1 , respectively. Similarly, HST images do not reveal 
any differences in the host galaxy magnitudes of these blazar 
subtypes either (Urry et al. 2000). 

The HST imaging survey of BL Lac objects by Scarpa et al. 
(2000) and Urry et al. (2000) demonstrated that the host galax- 
ies appear to be morphologically normal ellipticals. The mor- 
phological and dynamical properties of elliptical galaxies are 
linked via the fundamental plane (FP) relations (Djorgovski & 
Davis 1987; Dressier et al. 1987). Therefore, by combining the 
velocity dispersions with magnitudes and effective radii, we can 
test whether the BL Lac host galaxies fall on the same FP as 
normal ellipticals. For consistency with Urry et al. (2000) and 
Bettoni et al. (2001), we use H Q = 50 km s" 1 Mpc" 1 . 

One source of uncertainty in such a comparison comes from 
the fact that the effective radii measured from ground-based im- 
ages and from the HST snapshot survey of Scarpa et al. (2000) 
and Urry et al. (2000) are in severe disagreement for some of 
the galaxies in our sample. Table 4 lists the effective radii from 
the HST survey, and from ground-based imaging from various 
sources in the literature. The ground-based measurements of 
r e are nearly always larger than the HST measurements. In the 
worst case (3C 371), the disagreement is a factor of 5. 

The main advantage of the HST images is the narrow PSF 
core which makes it possible to determine the point-source 
flux accurately, but the HST measurements are hindered by 
the small field of view of the WFPC2 Planetary Camera (PC) 
detector (35" square), and brief exposures of 200 to 1300 s. 
High-quality ground based imaging (e.g., Nilsson et al. 1999) 
demonstrates that the host galaxies of the nearest objects are in- 
deed large enough to completely fill the PC chip. As a result, 
the sky background level determined from the outer portions of 
the detector will be contaminated by galaxy light, and measure- 
ments of the galaxy radial profile after subtraction of an artifi- 
cially high sky background would lead to an underestimate of 
r e . This systematic effect could be largely responsible for the 
differences between the ground-based and HST measurements. 
The total host galaxy magnitudes are in better agreement, with 
differences of typically less than a few tenths of a magnitude. 
The problem of determining the host galaxy parameters for BL 
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Lac objects from shallow HST images is somewhat analogous 
to the difficulty of detecting QSO host galaxies in short WFPC2 
exposures, as discussed by McLeod & Rieke (1995). Neverthe- 
less, it is clear that the ground-based measurements of r e are 
quite uncertain as well. The case of Mrk 501 illustrates the dif- 
ficulty of measuring r e from the ground. Even though this is one 
of the nearest and best-resolved BL Lac objects, ground-based 
measurements of r e range from 9" to 20". 

Figure 4 plots the FP for the BL Lac objects using the param- 
eterization given by Bettoni et al. (2001) for low-redshift radio 
galaxies. Given the uncertainty in the measurements of r e , and 
the consequent uncertainty in (fi e ) (the average /?-band surface 
brightness enclosed within r e ), we plot the BL Lac objects us- 
ing both the ground-based and HST measurements. The host 
galaxy magnitudes have been corrected for Galactic extinction 
and fc-corrected as described by Urry et al. (2000). The ground- 
based data included in the plot are from Abraham, McHardy, 
& Crawford (1991), Wurtz, Stocke, & Yee (1996), and Nilsson 
et al. (1999). For comparison, the radio galaxy sample of Bet- 
toni et al. (2001) and the Coma cluster ellipticals observed by 
J0rgensen, Franx, & Kjaergaard (1995a,b) are also shown. 

The overall agreement is very reasonable, considering the 
relatively large uncertainties in the measurements of a and r e 
for the BL Lac objects. Changes in r e at fixed total magnitude 
result in displacement along a vector that is nearly parallel to 
the FP sequence; as a result both the ground-based and HST- 
based measurements appear consistent with the FP. In either 
case, most of the BL Lac objects fall within the scatter of the 
FP. The largest outlier is 3C 371 when plotted with the effective 
radius of 2.9 kpc from Urry et al. (2000); it deviates from the FP 
sequence by ~ 0.3 dex. However, ground-based measurements 
indicate an effective radius that is larger by a factor of five, and 
with the ground-based value of r e this galaxy falls nicely on the 
FP. Thus, despite the uncertainty in the effective radii, the host 
galaxies do appear to be normal elliptical galaxies. 

The contrast between the ground-based and ffi>r-based mea- 
surements of r e and (fj, e ) can be seen more clearly by plotting 
projections of the FP onto the r e and (/z e ) axes, as shown in 
Figure 5. With the HST measurements, the BL Lac objects 
have effective radii comparable to the smallest radio galaxies, 
and most of the BL Lac hosts have (/j, e ) brighter than 19 mag 
arcsec" 2 , in a region where very few radio galaxy hosts lie. The 
ground-based measurements place the BL Lac hosts in a region 
that overlaps much more closely with the radio galaxies, bet- 
ter matching expectations from unified models. However, the 
BL Lac objects are still displaced toward smaller radii than the 
largest radio galaxies; this is also apparent in Figure 4, where 
the BL Lac objects do not extend as far up the FP sequence as 
the radio galaxies. Similarly, Urry et al. (2000) noted that no 
BL Lac hosts are found to be as luminous as the most luminous 
hosts of FR I radio galaxies. This offset could be a consequence 
of selection effects, or might indicate an environmental effect 
disfavoring the formation of blazar jets in rich clusters (Urry et 
al. 2000). 

The results presented here, and by FKT, strongly disagree 
with several recent estimates of M BH in blazars based on other 
methods. Using variability timescales to estimate the size of the 
emission region, Xie et al. ( 1 998) and Fan, Xie, & Bacon ( 1 999) 
derive estimates of the black hole masses for several blazars that 
are in the range ~ 10 6 - 10 8 M Q , much smaller than the masses 
obtained from the M BH - cr relation. Wang, Xue, & Wang (2002) 
have proposed a method to derive M BH from the peak lumi- 



nosity and frequency of the blazar emission. They conclude 
that HBL and LBL blazars have a bimodal distribution of M B h, 
with masses of 10 5 - 10 8 M for HBL and 10 9 - 10 11 M Q for 
LBL objects. The stellar velocity dispersions imply black hole 
masses that are are completely inconsistent with this conclu- 
sion. Models that interpret periodicities in blazar light curves 
in terms of orbital motion of black hole binary systems (Rieger 
& Mannheim 2000; De Paolis, Ingrosso, & Nucita 2002) also 
give much smaller masses than those derived from the Mbh - f 
relation. 

Xie et al. (2002) propose a new method to estimate black 
hole masses in BL Lac objects, based on the assumption that 
the observed variability timescale is determined by the orbital 
timescale near the innermost stable orbit around the black hole. 
They find masses ranging from 10 7 4 to 10 9 2 for 13 objects, in 
the same range as the masses we derive. Xie et al. show that 
their M BH estimates are strongly correlated with the beamed 
nonthermal luminosity but not with the host galaxy luminosity. 
However, the opposite would be expected if BL Lac host galax- 
ies follow the M B H-ibui g e correlation of normal galaxies. Fur- 
thermore, Urry et al. (2000) show that there is no correlation 
between the nuclear and host galaxy luminosities for BL Lac 
objects, implying that the nonthermal luminosity should not be 
tightly correlated with M BH if M Bli scales with host galaxy lu- 
minosity. Only two objects are in common between our sample 
and Xie et al., so it is not possible to perform a detailed com- 
parison of the mass estimates, but we note that their upper limit 
to the black hole mass in 3C 371 is an order of magnitude lower 
than the mass we derive from the M BH - a correlation. 

Wu, Liu, & Zhang (2002) and Woo & Urry (2002) have re- 
cently published estimates of the BH masses in BL Lac objects 
by using the HST imaging survey results to obtain approximate 
velocity dispersions via the FP relation, and then applying the 
M BH -o- relation to estimate M BH . The advantage of this method 
is that it can be easily applied to a large number of objects, as 
there are over 60 BL Lac objects with morphological parame- 
ters measured from the HST survey. However, there are some 
potential disadvantages as well, since the a estimates are lim- 
ited by the accuracy of the host galaxy decompositions and by 
the intrinsic scatter in the FP. In both cases, the FP-derived esti- 
mates of a are systematically larger than our measured values. 
For objects in common with our sample, the Wu et al. estimates 
of a are higher than our values by 30% on average, and the Woo 
& Urry estimates are higher by 74%. The most severe disagree- 
ment is for 3C 371; for this object Wu et al. predict a = 524 km 
s" 1 and Woo & Urry predict 618 km s" 1 , in comparison with our 
measured value of 249 ± 25 km s" 1 . This illustrates the need for 
accurate values of r e and (^ e ) when applying this method. 

Determining the Eddington ratio L/Lndd in BL Lac objects 
can give important clues to the physical structure of their ac- 
cretion flows. It would be interesting to search for systematic 
differences between LBL and HBL objects; this could yield 
new constraints for models of the spectral energy distributions 
(e.g., Ghisellini et al. 1998; Bottcher & Dermer 2002). Un- 
fortunately, it is particularly difficult to measure the bolometric 
luminosity of the unbeamed AGN core in BL Lac objects that 
are dominated by the beamed nonthermal emission from the jet. 
One alternative is to use the emission lines as a measure of the 
ionizing luminosity, because the emission-line flux is the only 
unbeamed component of the AGN that can be directly observed. 
This method is not ideal since there is no single emission line 
that can be observed in all objects; in practice the only viable 
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alternative is to use the lines that are detected and perform a 
rough conversion either to a total broad-line region luminosity 
(Celotti, Padovani, & Ghisellini 1997) or to a single line such as 
H/3 (Cao 2002). For our sample, broad-line fluxes are available 
only for a few objects, while a few others have [O III] emission 
detected. Given this small number of objects and the large un- 
certainties involved in comparing the luminosities of different 
emission lines, we are unable to perform any meaningful com- 
parison of the Eddington ratios of the HBL and LBL objects in 
our sample. This would be a useful objective for future work if 
stellar velocity dispersions and emission-line luminosities can 
be measured for larger samples of BL Lac objects. 

6. CONCLUSIONS 

Our conclusions are summarized as follows: 

1 . Stellar velocity dispersions in low-redshift BL Lac objects 
are in the range ~ 170-370 km s -1 . Using the Tremaine et al. 
(2002) fit to the Mbh - & relation, the corresponding black hole 
masses are M B h ~ 10 7 9 to 10 92 M Q . 

2. Measurement of velocity dispersions in BL Lac objects re- 
quires particular care in selecting spectral regions that are sensi- 
tive to a but are not severely affected by interstellar absorption 
lines, telluric absorption bands, or emission lines. Due to the 
dilution by nonthermal emission, high S/N (typically > 100 per 
pixel in the extracted spectrum) is required for accurate mea- 
surements with a direct template-fitting method. 

3. The distribution of velocity dispersions in BL Lac objects 
appears superficially similar to that of the radio galaxy sam- 
ple of Bettoni et al. (2001), indicating a similar distribution of 
black hole masses for the two samples. However, neither sam- 
ple is statistically well defined or complete in any sense and we 
caution against drawing conclusions about radio-loud AGNs in 
general from this limited comparison. 

4. There does not appear to be any systematic difference be- 
tween the black hole masses in HBL and LBL objects. 



5. The host galaxies of BL Lac objects lie on the fundamen- 
tal plane of nearby elliptical galaxies. We do not find any BL 
Lac objects near the top of the fundamental plane sequence in 
the region occupied by the most luminous hosts of FR I radio 
galaxies; this may be a consequence of selection effects in the 
identification of BL Lac objects. 

6. For nearby, well-resolved BL Lac objects, the host galaxy 
effective radii measured in the HST WFPC2 snapshot survey 
of Scarpa et al. (2000) and Urry et al. (2000) are systematically 
smaller than effective radii measured from deeper ground-based 
images. We suggest that the measurements of r e from the HST 
snapshot survey images may be biased toward underestimates 
of the true radii, because they were derived from shallow ex- 
posures that did not detect the faint outer envelopes of the host 
galaxies that are seen in deep ground-based images. 
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Table 1 
Observing Runs 



Observing 
Run 


UT Date 


Telescope 


Slit 

(") 


Grating 
(grooves mm" 1 ) 


A 

(A) 


Plate Scale 
(A pixel" 1 ) 


cr(instrumental) 
(kms" 1 ) 3 


1 


2001 Jun 24 


P200 


2 


Blue: 600 


4210-5950 


1.72 


115 










Red: 1200 


8400-9070 


0.63 


25 


2 


2001 Nov 9 


du Pont 


1 


832 


8240-9800 


0.72 


25 


3 


2002 Jan 22 


P200 


2 


Blue: 600 


4600-6350 


1.72 


115 










Red: 1200 


8400-9070 


0.63 


25 


4 


2002 Jan 24-26 


P200 


2 


Blue: 1200 


4940-5820 


0.88 


60 










Red: 1200 


8400-9070 


0.63 


25 


5 


2002 Feb 19-20 


Baade 


1 


600 


3600-6700 


1.56 


91 



a Instrumental dispersion for a source uniformly filling the slit, measured from the widths of comparison lamp lines at 5500 A or 8600 



A. 



Table 2 
Observations and Results 



Object 


Type 


z 


Run 


Exposure 


Extraction 


Seeing 


S/N 


a 


log(M BH /M Q ) 










(s) 


(") 


(") 




(kms" 1 ) 




Mrk 180 


H 


0.045 


4 


5400 


5.0 


1.5-2.0 


130 


209 ± 1 1 


8.20 ±0.11 


Mrk421 


H 


0.030 


3 


10800 


5.0 


2 


310 


219 ±11 


8.28 ±0.11 


Mrk 501 


H 


0.034 


1 


12600 


3.7 


1.5 


305 


372 ±18 


9.21 ±0.13 


AP Lib 


H 


0.049 


4 


4850 


5.0 


2 


160 


196±21 


8.09 ±0.20 


1 Zw 187 


H 


0.055 


4 


5400 


5.0 


1.5-2.0 


100 


171 ±12 


7.86±0.14 


3C 371 


L 


0.051 


4 


3600 


5.0 


3-5 


110 


249 ± 25 


8.51 ±0.19 


0521-365 


L 


0.055 


5 


1800 


2.5 


0.7-1.0 


165 


269 ± 29 


8.65 ±0.20 


0548-322 


H 


0.069 


5 


3600 


2.5 


0.7-1.0 


130 


202 ±24 


8.15 ±0.22 


0706+591 


H 


0.125 


3 


10800 


5.0 


2 


160 


216 ±23 


8.26 ±0.20 


2201+044 


L 


0.027 


2 


7200 


4.2 


1.0-1.5 


50 


197 ±8 


8.10±0.09 


2344+514 


H 


0.044 


3 


5400 


5.0 


2 


170 


294 ± 24 


8.80±0.16 



Note. — Blazar type is H = high-frequency peaked, or L = low-frequency peaked. Classifications are from the catalog of 
Scarpa et al. 2000, and redshifts are from NED. S/N is the average signal-to-noise ratio per pixel in the extracted spectrum, over 
the region used to measure the velocity dispersion. For Mrk 501, the S/N value is listed for the red side exposure from Paper I. 
Black hole masses are calculated using the fit to the M Bll - a relation by Tremaine et al. (2002). 
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Table 3 



Nearby Galaxies 



Galaxy 




a (km s" 






Red 


Blue 


McElroy (1995) 


NGC 2841 


222 ±4 


229 ±5 


232 ±18 


NGC 4278 


261 ±8 


251 ±13 


250 ±15 


NGC 4321 


92±4 


101 ±12 


94 ±12 


NGC 4374 


308 ±7 


302 ±7 


296 ±13 


NGC 4579 


165 ±4 


166 ±8 


189±4 


NGC 4736 


112±3 


109 ±4 


136 ± 12 



Note. — Data are from the June 2001 and January 2002 Palomar runs. 
Comparison results from the literature are taken from the compilation 
by McElroy (1995). The uncertainty ranges quoted by McElroy repre- 
sent the scatter among measurements taken from the literature for each 
galaxy. 



Table 4 



Effective Radii of BL Lac Object Host Galaxies 



Galaxy 


r e (HST) 


r e (ground) and Reference 




(arcsec) 


(arcsec) 


Mrk 180 


3.10±0.02 


7.1 (A91) 


Mrk421 


3.95 ±0.05 


1 1 (N99) 


Mrk 501 


4.44 ±0.03 


9.3 (A91), 20.0 (S93), 9.0 (W96), 13 (N99), 17.2 (P02), 15 (FKT) 


APLib 


3.70±0.10 


5.7 (A91), 8.7 (S93),6.7 (P02) 


lZw 187 


3.15 ±0.05 


3.9 (W96) 


3C 371 


2.10±0.10 


10.5 (W96), 14.9 (S93), 10.8 (P02) 


0521-365 


2.80 ±0.07 


3.4 (W96) 


0548-322 


7.05 ±0.15 


7.7 (W96) 


0706+591 


3.05 ±0.07 




2201+044 


6.78 ±0.08 


6.4 (F96), 7.6 (W96) 


2344+514 


5.93 ±0.02 


8.5 (N99), 5.8 (F99) 



Note. — HST measurements are from Urry et al. (2000), except for Mrk 501 which is described in 
§3.2. Ground-based measurements are taken from the following sources. A91: Abraham, McHardy, 
& Crawford (1991). S93: Stickel, Fried, & Kiihr (1993). W96: Wurtz, Stocke, & Yee (1996). 
F99: Falomo & Kotilainen (1999). N99: Nilsson et al. (1999). FKT: Falomo, Kotilainen, & Treves 
(2002). P02: Purismo et al. (2002). 
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FIG. 1 . — Best fits of broadened stellar templates to the BL Lac spectra. The thick curves show the broadened templates over the spectral region used to compute 
X 2 . In regions excluded from the fit, the broadened template spectra are shown with a dotted line. Galactic Na D absorption in the spectrum of 0706+591 is labelled. 
Mrk 501 is shown in Paper I. 
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FIG. 2. — Results of Monte Carlo simulations for measurement of the velocity dispersion from the diluted, noise-added spectrum of NGC 2841. A power-law 
continuum was added to contribute 2/3 of the total flux at 5500 A, and Poisson noise was added to yield a given S/N per pixel. Upper panel: The weighted average 
velocity dispersion measured from 100 realizations of the noise-added spectrum at each S/N value, with 10 template stars used for each measurement. Without 
dilution or added noise, the measured velocity dispersion is 229 ± 6 km s -1 . Lower panel: The standard deviation of velocity dispersions measured from the 100 
trials at each S/N value. 
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FIG. 3. — Upper panel: Histogram of stellar velocity dispersions in BL Lac objects. Cross-hatched points indicate the low-frequency-peaked objects. Lower 
panel: Low-redshift radio galaxies from the extended sample of Bettoni et al. (2001), which compiles data from several radio galaxy samples in the literature. The 
black hole mass scale at the top of the figure is based on the fit to the Mbh _ a relation by Tremaine et al. (2002). 
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FIG. 4. — The fundamental plane using the parameterization given by Bettoni et al. (2001). Filled squares: BL Lac objects, using effective radii and magnitudes 
from the HST survey of Urry et al. (2000). Open squares: BL Lac objects using effective radii and magnitudes from ground-based data. No ground-based 
measurments are available for object 0706+591. Crosses: Radio galaxies from Bettoni et al. (2001). Small circles: Coma cluster ellipticals from J0rgensen, Franx, 
& Kja;rgaard (1995a,b). 
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FIG. 5. — Projections of the fundamental plane onto the logc vs. logr c and (fi e ) planes. Filled squares: BL Lac objects, using effective radii and magnitudes 
from the HST survey of Urry et al. (2000). Open squares: BL Lac objects using effective radii and magnitudes from ground-based data. Crosses: Radio galaxies 
from Bettoni et al. (2001). 



